Anisotropic phonon DOS: the application of Rietveld 
and Mossbauer texture analysis in aligned powders 
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While synthesizing the single crystals of novel materials is not always feasible, orienting the 
layered polycrystals becomes an attractive method in the studies of angular dependencies of inelastic 
scattering of x-rays or neutrons. Putting in use the Rietveld analysis of layered structures in novel 
manganites and cuprates we develop the studies of their anisotropic properties with oriented powders 
instead of single crystals. Densities of phonon states (DOS) and atomic thermal displacememts 
(ATD) are anisotropic in the A-site ordered manganites LnBaMn20 y of both y = 5 and y = 6 
series (Ln=Y, La, Sm, Gd). We establish the angular dependence of DOS on textures of arbitrary 
strengths, link the textures observed by x-ray and 7— ray techniques, and solve the problem of 
disentanglement of Goldanskii-Karyagin effect (GKE) and texture in Mossbauer spectra. 

PACS numbers: 63.22.Np;76.80.+y 
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Two states of matter traditionally used in structure 
analysis are single crystals and polycrystalline powder. 
When synthesizing the single crystal is a stiff task the 
random polycrystalline materials or even polycrystals ori- 
ented on a surface conformably to their persistent crys- 
tal habitus can be employed. Via the method devel- 
oped initially by Rietveld[l| the refinement of texture, or 
preferred orientation, is a conventional procedure along- 
side with the refinement of atomic structure parameters. 
Randomization of the powders favors the perfect refine- 
ments of atomic parameters, therefore, any residual tex- 
ture is usually not in line with the best sample prepara- 
tion for x-ray profile analysis. There appears, however, a 
class of problems, in which the well-oriented powders of 
platy or acicular crystallites can replace the unavailable 
single crystals. To them belong the studies in anisotropic 
properties of materials, in electric, magnetic properties, 
and in lattice dynamics. The texture descriptions and de- 
termination of the preferred orientation parameters are 
thus becoming the issues of self-sustained interest. 

Recently, orienting powders to make the samples with 
varied degree of crystallite alignment was suggested Q to 
be useful in vibrational spectroscopy of anisotropic mate- 
rials. When the polycrystalline material is a ferromagnet 
or a superconductor it can be thoroughly subjected to 
texturing in an external magnetic field. 

The vibrational properties in novel manganites and 
cuprates are anisotropic owing to their layered structure. 
Substitution of Fe into 3d Jahn- Teller cation site perturbs 
the electronic system at small doping rate, and modifies 
the existing charge and/or orbital order. However, due to 
strong electron-phonon interactions the system keeps its 
attractiveness for high-resolution spectroscopies based on 
nuclear 7— resonance. The Fe-doped samples allowed us 




FIG. 1: Structure of oxygen-deficient layered oxides, e.g., 
manganites LnBaM^Os (a); geometry of nuclear inelastic 
scattering of synchrotron radiation (b); phonon density of 
states of Fe in unoriented powder of TbBaFe2 0s (c). 



to probe the anisotropy of the dopant vibrations due to 
GKE in 57 Fe Mossbauer spectroscopy 0. The same dop- 
ing procedure is extremely valuable in the synchrotron 
radiation vibrational spectroscopy using the nuclear in- 
elastic scattering (NIS) Q. In this work, we establish the 
relationships between the parameters of preferred orien- 
tation, the anisotropic phonon DOS, and the Mossbauer 
line intensities. Using these relationships we propose the 
novel technique of NIS on the oriented powder samples. 

Previous works dealt with the texture manifesta- 
tions in Mossbauer spectra [1, 0, 0]. For a unpolarized 
Mossbauer source there occur two parameters of texture, 
which can be determined from spectra. These parameters 
define the so-called "minimum texture function" (MTF). 
Our approach is to employ the March-Dollase[8( func- 
tion (MDF) implemented in programs for diffraction pro- 
file analysis [il lid]. We establish the relation between 
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Rietveld and Mossbauer textures in terms of MDF and 
MTF, using the example of quadrupole doublet spectra 
for the 57 Fe nuclei occupying 2% of Mn sites in param- 
agnetic state of novel manganites LnBaMn 2 05 [ll| and 
LnBaMn 2 6 [H H . The symmetries are tetragonal, 
except a monoclinic member YBaM^Og, and measure- 
ments were made with azimuthal rotation. Our texture 
function T{9) was uniaxial, (^—invariant (Fig.l). 

In the proposed method, the DOS derivation is based 
on two NIS patterns and a Rietveld pattern. Owing to re- 
lationship between MTF and MDF, the Mossbauer spec- 
tra and Rietveld analysis become mutually complemen- 
tary techniques. Indeed, the transmission (Mossbauer) 
and reflection (Bragg or NIS) data may diverge if 
there occur some in-depth variations of texture. The 
7— resonance wavelength (0.86 A) is close to x-ray Mo 
K a (0.7 A), but the main component of NIS radiation col- 
lected by the avalanche photodiode detector (APD) is Fe 
K a (1.94 A)[4], that is close to x-ray wavelength of Cu 
if Q (1.54 A). 

In the NIS spectrum of an anisotropic crystal, the 
phonon DOS is weighted by squared projection of the 
phonon polarization vectors to the wave vector of the 
x-ray quantum [l5|. Let the incident beam be launched 
under the angle "d with respect to preferred axis (z-axis) 
of a plate- or a needle-like crystallite. The projected DOS 
for this crystallite is: 



ar 1 / 2 /<& 



9eW = 9z{E) cos 2 + g x (E) sin 2 



(1) 



To introduce the averaging of the DOS over the ensemble 
of aligned crystallites we must integrate these two terms 
with the volume of crystallites D{d : <fi)dQ, whose z-axis 
lies within the cone shell clement dft: 
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FIG. 2: The r— dependence of the $— scaled coefficients of the 
Legendre series expansion of the MDF, Eq. [6] (a) and con- 
tour polar plot of the isotropic (a = 0) site $— scaled spectra 
asymmetry ARi/Q in function of the strength of preferred 
orientation r and incidence angle u) (b). 
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Prior measuring the NIS spectra Mossbauer spec- 
troscopy and Rietveld analysis provide the results in 
terms of MTF and MDF, respectively. Due to the ax- 
ial symmetry, the texture function T{6) would have only 
even terms in the Legendre expansion series. 
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(g(E)) = 9x (E) 



Ag zx (E) J D(i?,</>) cos 2 VdQ 



(2) 



Here Ag zx (E) = g z (E) — g x {E) and the orientation dis- 
tribution function (ODF) is normalized to unity. The 
polar ODF T(9) is to replace D{&, 4>) via the coordi- 
nate transform from the frame of the beam to the frame 
of the rotation stage. The ratio of angular elements 
cK^beam/rf^stagc is sin i)di)d(f>/ sin9d9dip and the Jaco- 
bian of this transform is sin ■d/sm.9. Using cos$ = 
cos 9 cos lo — sin 9 sin lo cos tp we obtain for the uniaxial 
symmetry 



The observable quantities, (g(E)), and asymmetries in 
Mossbauer spectra (R or ARi, as denoted below) are 
fully determined by the integrals of T(9), i.e., Q and 
(cos 2 0) = 1 — (sin 2 d), therefore, the function T[ff) con- 
sistent with experiment is not unique. The simplest func- 
tion T(9) conforming the Eq. |gj|, is the so-called MTF 
Q, composed of a linear combination of Pq = 1 and 
P2(cos6*) terms, i.e. .the par abolic function of cos 9. 

Rietveld analysis [H 0, H3| specifies the preferred ori- 
entation with two fit parameters r and $ of the MDF 
1: 



M(9,<f>,r) = 1 - $ 1 - ( 
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cos 2 uj + - (sin 2 9) sin 2 u (3) 



Powder averaging (sin 2 (#)) = J sin 2 (9)T (9) sin 9 d9 can 
be expressed via Legendre function of the first kind 
P 2 (x) = i(3x 2 - 1): 



The angle 9 describes the misfit of the crystallite ori- 
entation with respect to the axis of the uniaxial tex- 
ture (Fig.l, b), and <& is the fraction of the oriented 
phase. The March variable r expresses the strength of 
the preferred orientation. In a sample made by press- 
ing a layer of plates (needles) of the initial thickness do 
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(d) down to the thickness d (do), the degree of com- 
pression is r = d/dn(l6|. The function M(9, $,r) con- 
serves the scattering matter, therefore, it has a clearer 
physical sense than T(9) originally used by Rietveld 
T(9) = exp(— r-6> 2 )[l|]. MDF is thus a true normalized an- 
gular distribution, equally good to the textures made of 
crystallites with either platy or acicular habitus [U foL [lot • 
Ericsson and Wappling have previously studied in 
chlorite [ll] the effect of texturing induced by compres- 
sion of various degrees = di/do on the Mossbauer line 
area asymmetry. They have pointed out that their model 
"represent an extreme case" and that less randomness 
were observed in chlorite samples than predicted by their 
one-parametric compression model. Nagy [17] has intro- 
duced a second parameter taking into account the par- 
ticle shape and assumed that the texturing behavior of 
the particles under compression is shape-dependent. For 
the platy chlorite flakes Nagy obtained the aspect ratio 
(cylindrical h/D) as large as 0.6 (l7j . The March-Dollase 
function implemented in the FULLPROF program [l0| is 
also two-parametric, however, the meaning of the second 
parameter $ is unrelated to the method of particle align- 
ment. The unaligned random phase is just added with 
the fraction of 1— $. We show below that the Mossbauer 
line area asymmetry depends for small asymmetries only 
on the product <frr. Difference of physical meanings be- 
tween $ and r manifests itself only in strongly asymmet- 
ric spectra. 

Measurements with the rotated sample stage inclined 
to the incident beam (Fig.l, b) provide us with two com- 
ponents g x (E), g z (E) of uniaxial DOS in terms of the 
angle ui and the parameters r and <E>. Using MDF (EqJSJ), 
we obtain 

(sin 2 9) = V(r, $) = f M(6, $, r) sin 3 8d9 (7) 
Jo 

(8) 
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e(r) = \/r 2 — i — 1 
f3(r) = ln(2r 3 + 2y/r e - r 3 - 1) 



(9) 
(10) 
(11) 



(cos 2 0} = [1 - V(r, $)] cos 2 w + 



V{r,$). 2 



sin 2 w (12) 



via substitution (sin 2 #), (cos 2 9) into Eq.©. Both e(r) 
and (3(f) are imaginary for < r < 1, however, v(r) is 
real for < r < oo. The ranges < r < 1, and 1 < r < 
oo correspond to platy and stalky habits, respectively. 
Now from © we can express MTF through r and $: 

Co 5 _ _ . . /15 45 2 \ „. , 15 2 3 

t + 2 C2P2(a;) = U ~ - x ) V ^ $) + y x - 2 

(13) 

It is shown in Fig. 2 (a) that the Legendre series coeffi- 
cients of the MDF are (all, except Co = 1) returning to 
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FIG. 3: Vibrational anisotropy a = k 2 ((z 2 ) — (x 2 )) in func- 
tion of three variables a(r, &,ARi) , Eq. [22], is shown as 
three slices (contour plots) of a(r, <&). In each ARi —matched 
slice, the oriented samples of manganites LnBaM^Os and 
LnBaMn2 05 are shown by Ln symbols and oxygen index as 
spots with coordinates r, $ obtained from Rietveld refinement 
of parameters of preferred orientation. 



zero at r = 1 and that the MDF is well approximated 
by \Cq + (5/2)6*2 P 2 (x) m some vicinity of random poly- 
crystal (r = 1), because C n ~ for n > 2 around r = 1. 

We are ready for determination of both g x {E) and 
g z (E). From Eq.© a couple of measurements of DOS 
gi(E) and gi(E) at the angles u>i and W2 leads to: 



Ag zx (E) = 



[l - § V(r, $)] (cos 2 wi - cos 2 w 2 ) 



(14) 



lTl — / rp\ AgiaCg) (cos 2 ^i +cos 2 cj 2 ) , , 

g x (E)=g(E) ^ 2 \ ( 15 ) 

2 (cos z wi — cos z W2J 

#(£) = gi(E)/2 + g 2 (E)/2 , A 5 i 2 (£) = 91(E) - g 2 (E). 

Now we turn to the asymmetry of Mossbauer spectra 
caused by texture with parameters r and First, the 
intensity ratio for a single crystallite is to be examined 
depending on the orientation of the wave vector of the 
incident x-ray quantum with respect to the axes of the 
electric field gradient (EFG) tensor at the site wherein 
the 57 Fe nucleus is located. The ^-dependent Clebsch- 
Gordan coefficients determine the doublet line intensity 
ratio [fjj]. We found that the quantity, which is pro- 
portional to the amount of the oriented phase <E> is the 
deviation of relative line area from 1/2: 
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Powder averaging resolves, using Eq.©, into 
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FIG. 4: Mossbauer spectra in oriented samples of 57 Fe-doped 
manganites. Two panels present the oxygen-poor and oxygen- 
rich series ("O5" and "OV')- The spectra are fitted either 
with one or with two asymmetric doublets. Meltdowns of 
Mn 2+ /Mn 3+ and Mn 3+ /Mn 4+ charge order in 05 and 06 
series, respectively, are seen as single-site doublets. The 
occurrence of two sites for Fe in La and Sm members of 
the LnBaMn2 05 series is interpreted as remainder of un- 
molten charge order between Mn 2+ and Mn 3+ . The asym- 
metries A7?i were fixed to be a single parameter of both 
doublets. This is because each crystallite contain both these 
sites and their V zz direction and sign coincide according to 
structurelllll. 
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Eq. (fl"7f remains true in $— scaled form, i.e. with ARi /$ 
in place of ARi, and v(r) in place of V(r, $). From Fig. 
2(b) starting with ARi /$ one finds the MDF variable r 
at any angle co except magic angle u>m — 54.7°. 

When the measurement is done at u> = ujm the texture 
does not manifest itself in the spectra. In a simple normal 
incidence transmission Mossbauer spectrum (u = 0), the 
expression for intensity ratio was formulated [5(, 



vector times the difference of mean-square vibrational 
displacements along V zz , and in perpendicular direction, 
a = k 2 ({z 2 ) - (x 2 )), k 2 = 53.35A" 2 . Here again, it is 
more convenient to work with ARi = 1(1 — R)/(l +R). 
The replacement of M(0,$,r) with the MTF, Eq. [13l 
makes our expression for Ai? i integrable and expressible 
via two simple sigmoid functions erg (a) and <Ti(a): 



ARi = 



_ <T (a) + $U{r) [fa x (a) - 5<r (a) - £ 



l-20$U(r)a (a) ^ 
with U(r) = 1 - 3v(r)/2, <n(a) = (1 + 3/2a)cr (a), and 
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Here K{a) = xFx(h, §, —a) is the Kummer confluent hy- 
pergeometric function. The range of variation for both 
sigmoid functions 00(a) and <ri(a) is between —1/4 and 
+1/8; for unoriented powder r = 1, v(l) = 0, and 
Aiii(a) = <to (a) is the exact solution. The accuracy 
of the solution (fT9")) is better than 1% only in the narrow 
range of slight textures (0.9 ~ ?' ~ 1.1 and < <& < 1), 
in which the MTF is a perfect approximation for MDF 
(Eq. [6]). In the same range, the Ai?i(a) differs from 
(Tq (a) merely by shift of sigmoid flcxpoint: 



ARi(r, a) = a (a - a (r, $)) 



(21) 



The approximation (|2lj) is as good as (fT9j) for ao = 
4.53(1 — $r). The parameters r and $ enter to ao in 
equivalent form, however, become nonequivalcnt if the 
o — approximation (|21|) is extended to a broader range 
of r. We extend it by taking into account the nar- 
rowing sigmoid for stronger textures: Ai?i (r, $, a) = 
ao(A(r, $)a — B(r, $)). The accuracy of this approxi- 
mation is better than 1% in the broad range of textures, 
0.3 is r ;$ 2. The r and <I> dependencies of A and B can be 
well fitted withA(r, $) = l+$Ai(r)+$ 2 A 2 (r) + $ 3 A 3 (r) 
and B(r,U) = <S>B 1 (r) + $ 2 B 2 (r) + <P 3 B 3 (r). The polyno- 
mials Ai{r) and Bi(r) could be expressed as dot products 
of coefficients vectors (ao, ax, a?, as, 04), with the vector 
(1, r, r 2 , r 3 , r 4 ) *18|. In practice, it is crucial to find the 
vibrational anisotropy a starting from ARi, therefore, 
a is given using the function inverse of 00: 



J V2 M(9, $, r) (1 + cos 2 6>)e~ Q 
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however, the solutions were yet found either for texture 
effects, or for GKE, separately only. In the proposed syn- 
chrotron experiments on anisotropic powders, both are 
crucial, therefore, the combined effects of texture and 
GKE are of our interest. In (fT5)) a is the squared wave 



a(r,$, ARi) 



A(r, <&) A(r,$) 



(22) 



The function of 3 variables a(r, <&, Ai?i ) is visualized 
via several 2D maps (slices, Fig. 3). Total range, in which 
CTo 1 (ARi) is defined, is -0.25 < ARi < 0.125, and the 
slices are shown for ARi = -0.01, -0.025, and -0.04. 
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FIG. 5: The angular dependence of the doublet asymme- 
try (a) as given by the Eq. (17) and the Mossbauer spec- 
trum measured under the magic angle ujm = 54.7° between 
the direction from source to detector and the normal to the 
GdBaMn1.g6Feo.04O6 sample absorber. The distance between 
source and the sample center was L = 4 cm and the illu- 
minated area of the sample was 2x2 cm (d x d). The an- 
gles at the absorber surface were thus varied between u>\ = 
com - tan- 1 ( . d ? 8 . a, M ) and uj 2 = lj m + tarrV ^° s "m ), 
This range between uj\ =45° and U2 = 61° is shown in (a) 
by bold line. A small residual asymmetry might have been 
resulted from the averaging in this range. The fitted value of 
AR1/2 = 0.0045 is 5.5 times smaller than the normal-incidence 
value of AR4/2 in Fig. 4. 



its discrepancy reaches 10% already at ARi = —0.05. 
The parabolic model a^ 1 = 31.4A.Ri + 96.5Ai?f is bet- 
ter than 0.1% in the range -0.05 < ARi < 0.05 that 
covers the asymmetries observed in our spectra (Fig. 4). 

We have checked that the rotation of the sample 
around an axis perpendicular to the k —vector of the 
incident 7— quantum induces a difference of the dou- 
blet asymmetry. The same as in Fig. 4 absorber of 
GdBaMii! 96 Feo.o406 was installed under the magic an- 
gle u>m = 54.7° and the symmetric doublet was obtained. 
A very small residual asymmetry can be explained by the 
divergence of the incident beam (Fig. 5). 

Both oxygen-saturated and oxygen-depleted families of 
layered manganites show negative ARi, although their 
left-to-right area ratio changes the sign between "O5" 
and "Oq" series. This is because the main axis of EFG 
(z) is perpendicular to layers and V zz > in the pyramid 



FeOs, but V zz < in the FeOe octahedron compressed 
along z— axis. Since the sign of V zz changes from "O5" to 
"0 6 " series the doublet lines swap their positions. ARi 
remains negative according to its definition in Eq. (|16[) . 
Therefore, in both panels of Fig. 4, the stronger line is 
±3/2 — ► ±i, and the weaker line is ±i — > ±i. In this re- 
spect, the layered manganites are quite similar to layered 
cuprates, wherein the ionic point charge model prescribes 
the same orientation and sign of y^[19j]. 

The values of r and $ for our sa mples;20] refined from 
x-ray diffraction patterns[l(| are shown by spots in Fig. 3, 

at slices of the same AR 1 as shown in Fig. 4. The value of 

2 
a increases with increasing size of Ln within either "O5" 

or " Oe" series, however, the vibrational anisotropy is pos- 
itive in pyramid FeOs {(z 2 ) > (x 2 )), but negative in the 
compressed octahedron Fe06 ((z 2 ) < (x 2 )). This result is 
verifiable via refinement of the factors of anisotropic ther- 
mal displacements (ATD) of Mn from neutron diffraction 
data, however, till now such the data were refined with 
Bi SO — models 13, On the other hand, in a similar hi- 
layered structure of YBaFeCuOs, the ATD of Fe shaped 
in FeOs pyramid as prolate 'cigar' was found[2l|. Also, 
the oblate 'pancake' is not unexpected shape of ATD in 
the oblate octahedron of LnBaMn^Oe- 

In conclusion, we have shown how the information con- 
tained in the x-ray diffraction full-profile patterns of the 
aligned powders can be used for finding the direction- 
projected components of the phonon DOS. This infor- 
mation is sufficient, in principle, to solve the problem of 
determination of two DOS components from at least two 
experiments either performed on two samples with dif- 
ferent and non-zero degrees of alignment or conducted 
on one aligned-powder sample but under at least two dif- 
ferent angles u>. The same information is contained in 
assymetry of Mossbauer line intensities, however, here 
the texture effect is entangled with GKE. Therefore, the 
useful for DOS or for other properties information can 
be disentangled from Mossbauer spectra in only cases 
when either GKE or ODF are characterized by some sep- 
arate experiments. One possibility is to determine the 
GKE separately from the ATD refinement using neutron 
diffraction full-profile analysis superior in the accuracy 
of B an i SO compared to x-ray Rietveld analysis. Another 
possibility is to employ the texture March-Dollase pa- 
rameters. The latter possibility was realized in this work 
for two cases when the sign and orientation of uniaxial 
V zz is known, namely, in LnBaM^Os (V zz > 0) and 
in LnBaMn206 (V zz < 0). Our method is well appli- 
cable for these tetragonal structures, in which the axis 
of preferred orientation coincides with the direction per- 
pendicular to layers and with the principal axis of elec- 
tric field gradient. More complicated textured systems, 
such as non-uniaxial or magnetic textures were beyond 
the scope of the present work. It must be finally em- 
phasized that the scope of possible applications of the 
present work is limited by the close correspondence be- 
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tween the distribution of March-Dollase, which is theo- 
retically substantiated only in pressed samples [17|. and 
the real ODF T(0), which can be obtained experimentally 
with texture goniometers using several analytical meth- 
ods, best of which combine the ODF calculation with 
the Rietveld structure refinement [22j|. Our approach is 
also good for samples containing the randomly oriented 
phase, because we employed the extended (biparametric) 
MDF M(6, $, r). 

This work was supported by Asahi Glass Foundation 
and RFBR-JSPS (Grant 07-02-91201). 
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FIGURE CAPTIONS 

Fig.l. Structure of oxygen-deficient layered oxides, 
e.g., manganites LnBaMn 2 05 (a); geometry of nu- 
clear inelastic scattering of synchrotron radiation (b); 
phonon density of states of Fe in unoriented powder of 
TbBaFe 2 5 (c). 

Fig. 2. The r— dependence of the scaled coefficients 
of the Legendre series expansion of the MDF, Eq. [S] 

(a) and contour polar plot of the isotropic (a = 0) site 
$— scaled spectra asymmetry ARi/(5> in function of the 
strength of preferred orientation r and incidence angle uj 

(b) . 

Fig. 3. Vibrational anisotropy a = k 2 ((z 2 ) — (a; 2 )) in 
function of three variables a(r, $,Ai?i) , Eq. [22] , is 
shown as three slices (contour plots) of a(r, $). In each 
ARi —matched slice, the oriented samples of manganites 
LnBaMn205 and LnBaMn 2 Os are shown by Ln symbols 
and oxygen index as spots with coordinates r, <£> obtained 
from Rietveld refinement of parameters of preferred ori- 
entation. 

Fig. 4. Mossbauer spectra in oriented samples of 57 Fe- 
doped manganites. Two panels present the oxygen-poor 
and oxygen-rich series ("O5" and "Oq"). The spectra are 
fitted either with one or with two asymmetric doublets. 
Meltdowns of Mn 2+ /Mn 3+ and Mn 3+ /Mn 4+ charge order 
in 05 and 06 series, respectively, are seen as single-site 
doublets. The occurrence of two sites for Fe in La and 
Sm members of the LnBaMn 2 05 series is interpreted as 
remainder of unmolten charge order between Mn 2+ and 
Mn 3+ . The asymmetries ARi were fixed to be a single 
parameter of both doublets. This is because each crystal- 
lite contain both these sites and their V zz direction and 



sign coincide according to structure[ll|. 

Fig. 5. The angular dependence of the doublet asymme- 
try (a) as given by the Eq. (17) and the Mossbauer spec- 
trum measured under the magic angle com = 54.7° be- 
tween the direction from source to detector and the nor- 
mal to the GdBaMni.g 6 Fe .04O6 sample absorber. The 
distance between source and the sample center was L — 
4 cm and the illuminated area of the sample was 2x2 
cm (d x d). The angles at the absorber surface were 

dcos " M ^ and 



thus varied between w x = w M ~ tan" 1 ^ dcosuJ ™ ) 

- 1 iVi V L —d sin cjm ' 

lu 2 = lum + tan"^ r ^°.""„ ). This range between w x 
=45° and w 2 = 



^ L-\-d sin wm ' 

61° is shown in (a) by bold line. 



A 



small residual asymmetry might have been resulted from 
the averaging in this range. The fitted value of A~Ri/ 2 
= 0.0045 is 5.5 times smaller than the normal-incidence 
value of AR]y 2 in Fig. 4. 



